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1620 PROBIT ANALYSIS PROGRAM 



By Georgia Lee V. Meivin 
IBM Corporation 
2330 St. Paul Street 
Baltimore 18, Maryland 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce, 
ment occurs, users should order a complete new program 
from the Program Information Department. 



o 



'-.ITH-UP 



HE r'i'.OD : The method described in D. J. Finney, Probit Analysis , 

Cambridge University Press, is used in this program. The -ORTRAN 
notation follows as closely as possible that ussd in the text (e.g. 
FORTRAN "EiV is equivalent to Finney's "n" , SN".TX2 is equivalent to 
Dnwx- ; such changes should be fairly obvious). An initial esti- 
mation of a = 5, b = is used in the iteration. 



T 2 7t _ 



r z 2 
/ - JL 

' e 2 



dx 



is approximated by a modification or a formula appearing in Hastings, 
Approximations for Digital Computers . The approximation is : 

k .- 4 ■ -rr ± v ' 

i-1 



P = 



for t<0 



where u 



r 1 r~ ^ i 

1 - J e .72 /__/ a u for t > 

! 72T 

[i + q ItO" 1 



i=i 



aj * 0.22583685, a 2 = -0.25212867 

a - 1.2596951, a » -1.2878224 
3 * 

a 5 - 0.94064607, 9 - 0.23164189 

The probit slope, b, is iterated until the relative error 

(New b - old b/New b) < .0001 or until 13 iteration* are completed 
if the slope does not converge in fewer cycles. This accuracy may be 

altered by changing the constant in the statement ZPSLN = 1. E-4 or 

the maximum number of cycles may be changed by altering the statement 
DO 4 J - 1, 13. 



RESULTS : The program outputs include the probit. equation (Probit « 2 

-t b (log dose)), 
chi-square and the associated degrees of freedom, and the 
log of the expected dose, the expected dose, with 951 con- 
fidence limits on both, for a maximum of 10 levels (e.g. 



50 • 
data. 



ED ); 
98 



those levels chosen are a part of the input 



R£C0MPILATI0N : If it is desired to recompile the source program for use 
on a basic machine or to change any parameters the following restrictions 
must be observed. 

a. Two extra subroutines are required. ABS refers to an ab- 
solute value subroutine. HEAD refers to a special purpose 
subroutine for printing the long alphabetic heading line. 
Symbolic listings of both subroutines follow the listing 
for the main program. They must be added to the FORTRAN 
subroutine deck according to the directions supplied in 
the FORTRAN write-up. 



b. If the punched card output option is desired, FORTRAN must 
be modified to allow 80 columns of output. See 1620 
Technical Bulletin # 115 (Technical Digest No.: 1620- 
PRG-20.44). 

c. If the program is intended for a paper tape machine , output 
oust b* on the typewriter or the source program must be 
changed so that all "PUNCH" statements read "PUNCH TAPE" . 
All "READ" statements mist be changed to "READ TAPE" or 
"ACCEPT" . 



INPUT DATA : 

1. One card containing Aj- A $ in p. type notation. Thi» card 

is constant and must be used with every run. It la iaclmded 
with the prograa deck. 

2. One card containing M in cc 1-2. This is the number of P 
▼aloes for which ED computations are desired, if must be 
an integer < 10. 

3. M cards containing OP and OT (1 pair per card). "OP" 
is the "P" value on the printout and "OT" is its 
associated prebit. 

ce 1-3 OP 
cc 4-11 OT 



l®ta vtlati exist be in Floating-point, 

F-type aas&tsion. 
The input cards a!*a«ly described remain constant for 
any one run. Following these are the data cards for 
each assay. 

4. One card with title for assay in cc 1-20. 

This may be a maximum of 20 alphanumeric characters.. 



3 







5. 



One card containing N in cc 1-2. N is the number 
of dosage groups for the assay. N must be an integer 
£50. 



6. N cards containing the following information for each 
group: 

cc 1-3, D (The dosage level); 

cc 4-8, R (The no. of subjects reacting); 

cc 9-16, EN (the no. of subjects in the group). 

All three values oust be in floating- 
point, F-type notation. 

The data for any number of assays (4,5,6 above) say £•!!•«. 



OUTPUT: The format of the typed output is show with the sample problem. 
The typewriter is programmed to space to a new peg* before 
beginning the next assay. Card output ia an exact asallcate 
of typewriter output including blank cards within the body of 
the problem (but not the page spacing at torn ami) set that an 
80-00 listing on the 407 will produce tka amm ftfflffit as tie* 
typed output . 
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SAMPLE INPUT 



.22583680 .25212370 1.25969510 1.28782250 .94064610 

1. 2.6737 

16. 4.0055 

30. 4.4756 

50. 5.0000 

84. 5.99^5 

99. 7.3263 
PRCBIT EXAMPLE J 

4 

5. 1. 10. 

10. 3. 10. 

25. 8. 10. 

50. 10. 10. 

PRCBIT EXAMPLE 2 
5 

1. 1. 10. 

5. 1. 8. 

10. 1. 4. 

25. 10. 20. 

50. 35. 40. 
PRCBIT EXAMPLE 3 

5 

35. 9. 10. 

20. 14. 20. 

10. 4. 8. 

5. 6. 20. 

1. 2. 10. 
PRCBIT EXAMPLE 4 

5 

1.0 1. 10. 

1.5 4. 20. 

2.0 tO. 40. 

5.0 6. 10. 

10. 5. 5. 



a 



P3C3IT EXAMPLE 1 

PRCBIT = +1.119 +3.458LCG DCSE 

O.F. CHI -SQUARE 

+2 +.4792 



95 PER 


CENT CCNFI 


DENCE LIMITS 


(NCN-SIG. CHI- 


-SQ.) 






P 


LCG ED(P) 


LCWER 


UPPER 


ED(P) 


LCWER 


UPPER 


+1. 


+.4493431 


+.0893390 


+.8093472 


+2.8141 


♦1.2283 


♦6.4468 


+ 16. 


+.8344083 


+.6291276 


+1.0396891 


+6.8298 


+4.2572 


+10.9569 


+30. 


+.9703290 


+.8036265 


+1.1370315 


+9.3396 


+6.3624 


♦13.7098 


+50. 


+1.1219495 


+.9737812 


+1.2701178 


+13.2418 


♦9.4141 


+18.6259 


+84. 


+1.4094907 


+1.2086184 


+1.6103630 


+25.6738 


♦16.1665 


+40.7720 


+99. 


+1.7945559 


+1.4404166 


+2.1486952 


+62.3097 


♦27.5687 


+140.8300 
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PRCBIT EXAMPLE 2 



PRCS ! T EXAMPLE 3 



PRCBIT = +3.015 +1.679LCG DCSE PRCBIT = +3.767 +1.352LCG DCSE 

O.F. CHI -SQUARE D.F. CH I -SQUARE 

+3 +6.0823 +3 +2.1531 

95 PER CENT CONFIDENCE LIMITS (NCN-SIG. CHI-SQ.) 95 PER CENT CONFIDENCE LIMITS (NCN-SIG. CHI-SQ.) 

P LCG ED(P) LOWER UPPER ED(P) LOWER UPPER p LOG ED(P) LOWER UPPER ED(P) LOWER UPPER 

+1. -.2035168 -.9502609 +.5432272 +.6258 +.1121 +3.4932 + i. -.8086213 -1.8001502 +.1829074 +.1553 +.0158 +1.5237 

+16. +.5896205 +.1790390 +1.0002019 +3.8870 +1.5102 +10.0046 +]6. +.1761488 -.3140370 +.6663347 +1.5001 +.4852 +4.6380 

+30. +.8695828 +.5659155 +1.1732502 +7.4059 +3.6805 +14.9021 +30. +.5237539 +.1859431 +.8615647 +3.3400 +1.5344 +7.2705 

+50. +1.1818830 +.9707824 +1.3929836 +15.2013 +9.3493 +24.7163 +50. +.9115099 +.6724800 +1.1505399 +8.1566 +4.7041 +14.1429 

+84. +1.7741454 +1.5233253 +2.0249655 +59.4491 +33.3676 +105.9169 +8k. +1.6468711 +1.2090628 +2.0846794 +44.3477 +16.1831 +121.5288 

+99. +2.5672828 +2.0100974 +3.1244682 +369.2179 +102.3522 +1331.8895 +99. +2.6316413 +1.6992081 +3.5640745 +428.1947 +50.0274 +3665.0043 



o o _ ^ ; ; _ __ © 



PRCBIT EXAMPLE 4 



PRCBIT « +3.585 +2.658LCG DCSE 



D.F. 

+3 



CHI -SQUARE 
+1 .0789 



95 PER CENT CONFIDENCE LIMITS (HCN-SIG. CHI-SQ.) 



p 


LOG ED(P) 


LOWER 


UPPER 


ED(P) 


LOWER 


UPPER 


+1. 


-.3429391 


-.7032123 


+.0173341 


+ .4540 


+.1980 


+1.0407 


+16. 


+.1580888 


+.0127159 


+.3034616 


+1.4390 


+1.0297 


+2.0112 


+30. 


+.3349421 


+.22t11Q1 


+.4487741 


+2.1624 


+1.6638 


+2.8104 


+50. 


+ .5322233 


+.3824933 


+.6819532 


+3.4058 


+2.4126 


+4.8078 


+84. 


+.9063578 


+.5994673 


+1.2132482 


+8.0604 


+3.9761 


+16.3398 


+99. 


+1.4073858 


+.8590527 


+1.9557189 


+25.5497 


+7.2285 


+90.3064 



COmPl 

TECHNOLOGY 



a 

* 



OPERATING INSTRUCTIONS 

1. Set Console Switches 



cowsotE saiTca cm off 

1 Program pauses Program runs 
after each assay; continuously; no 
push start to con- last card test 
tinue 

2 Punched output Typed output 

3 Not used 

4 Not used 

2. For typed output, typewriter margins oust be set to 
allow for at least 80 characters. No tab stops are 
necessary. Set to double space. 



3. Set parity and I/O switches to stop, overflow switch 
to program. 



4. Clear Storage: 



a. Press reset 

b. Press insert 

c. Type 1600010 00000 

d. Press release and start 

e. When clear, press instant stop 

f . Press reset 



5. Place program deck in card reader, press load . 

6. When program is loaded, message "load data" will be typed. 
Place data package in card reader, press Reader Start and 
and Console Start . (If desired, the data may be placed 
directly behind the program when program is loaded. When 
"Load data" is typed, press Console Start only). 

7. If a number of assays are to be run switch 1 may be left 
off and turned on as last assay starts. If switch 1 is 
off when reader runs out of cards, machine will hang up 
with reader no feed light on. 

8. If it is desired to change values of OP between assays, 
the reader must be emptied of cards and the complete 
data package as described in Input data section placed in 
reader. Press reset , insert , type 49 03300, r elease and 
start. 



f? t rt o \ 

A, ft, ft, />, /U 



! I?.,,!': for £0 

\ <t,p o,-.-< or) 



Re « o 
t, -Me o-F 
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Y 
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32 




13 



rTfi 2 fin 


r 
L 


08300 








UOJi *+ 






1 00 






58462 


101 


58484 




58 if 96 




58616 


102 


(38644 

%/ V / w ~ "T 


2 


586^6 


103 


fiR7?n 


mikk 




58756 




58828 


104 


58860 




58944 3 


09040 




59064 




59088 




59100 




59124 




(59148 




59172 




591 96 




59220 




59244 




59256 




59328 




59424 




59492 


6 


59552 




59636 


23 


59792 


60 


59860 


8 


59896 


9 


59980 




Toil 2 




To 148 




T0220 




T0268 




T0376 




T0436 




T0520 


5 


T0556 




To 5 92 




T0628 




T0688 




T0748 




T0784 




T0844 




To 90 4 




T0928 




T0996 4 



ANALYSIS , ■ ■ . v ■ % , % 

ON X(50),P(50),OP(10),OY(10),EN(50) 
1.E-4 

A2, A3, A4, A5 
F11.8, F11.8, F11.8, F11.8J 



- 1 ,n 
12, 0P(l), 
(F3.0,F8. 



0Y(l) 
*0 



1,13 



PROS I T 
D I MENS 
EPSLN 

READ 100, A1, 
FORMAT (F9.8, 
READ 101, M 
FORMAT (12) 
DO 1 I - 1,M 
READ 102 
FORMAT 
READ 103, 
FORMAT (20H 
READ 101, N 
DO 3 I * 1,N 
READ 104, D, R, EN ( I ) 
FORMAT (F3.0,F5.0,F8.0) 
P(I)=R/EN(I) 

X( I ) = .43429^ * LOG(D) 
A . 5o 
B « 0. 
DO 4 J - 
SNW= 0. 
SNWX =0. 
SNWY=0. 
SNWXY«0. 
SNWX2=0. 
SNWY2«0. 
DO 5 I - 1,N 
T - A + B- * X(l ) 
TE a (T-5.)*(T-5.)*.5 
IF (TE-12.5) 6,5,5 
Z = .3989423 * EXP(-TE) 
U= 
PW 

IF (T-5.)9,8,8 

PW m 1 .-PW 

Y - T + (P(l )-PW)/Z 
W - Z*Z/(PW*(1.-PW)) 
SNW « SNW + W 
SNWX - SNWX + W*X(I ) 
SNWY » SNWY + W*Y 
SNWX2 * SNWX2 + W* X(l ) 
SNWY2 - SNWY2 + W*Y*Y 
SNWXY = SNWXY + W * X(l ) 
CONTINUE 
XBAR « SNWX/SNW 
YBAR * SNWY/SNW 
SXY = SNWXY - YBAR * SNWX 
SXX a SNWX 2 - XBAR * SNWX 
BNXT . SXY / SXX 
A « YBAR - BNXT * XBAR 
U a ABS( (BNXT-B )/BNXT) 
B a BNXT 

IF (U-EPSlN) 7,7,^ 
CONTINUE 



dPi; 



1. /(1.+. 23164189* ABS(T-5.)) 

= ((((A5*U-A4)*U+A3)*U-A2)*U+A1)*U*Z *1. 414214 



EN(I ) 



X(l) 
★ y 
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T1032 7 N = N-2 

T1068 CHI = SNWY2 - YBAR * SNWY - SXY * BNXT 

T1176 W = 1. / SXX 

Tl 212 SNW = 1. / SNW 

T1248 Y = 1./B 

Tl28i* IF (SENSE SWITCH 2) 70, 71 

T1304 70 PUNCH 103, 

T1316 PUNCH 105, A, B, . 

T1352 105 FORMAT ( /8HPR0BIT =F7.3,F7.3,8HLOG DOSE/) 

TH36 106 F0RMAT(25H D.F. CH I -SQUARE/ 1 3 , 1 3XF8. 4/ ) 

Tl56Jf PUNCH 106, N, CHI 

T1600 PUNCH 107, 

T1612 107 F0RMAT(48H95 PER CENT CONFIDENCE LIMITS (NON-SIG. CHI-SQ.)) 

Tl 732 GO TO 72 

T1740 71 TYPE 103, 

T1752 TYPE 105, A, B 

T1788 TYPE 106, N, CHI 

T1824 TYPE 107. 

T1836 72 U . HEAD(U) 

T1860 DO 12 I - 1,M 

T1872 D - (0Y(I ) - A ) * Y 

T1944 R - - XBAR 

T1980 SX - Y * SQRT(SNW+ R*R*W) * 1.96 

T2076 DU « D + SX 

T21 12 DL - D - SX 

T2148 DE m EXP (2.3025851*0) 

T2196 DUE » EXP (2.3025851*0U) 

T2244 OLE - EXP (2.3025851*DL) 

T2292 IF (SENSE SWITCH 2) 80, 81 

T2312 80 PUNCH 109, OP(I), D, DL, DU, DE, DLE, DUE 

T2432 109 F0RMAT(F4.0, F13.7, F13.7, F13.7, F12.4, F12.4, F12.4) 

T2484 GO TO 12 

12*92 81 TYPE 109. OP(l), D, DL, DU, DE, DLE, DUE 

1261 2 12 CONTINUE 

T2648 IF (SENSE SWITCH 1) 13, 35 

T2668 13 PAUSE 

T2680 35 IF (SENSE SWITCH 2) 37, 36 

T2700 36 H - 15 - M 

T2736 DO 30 I - 1, N 

T2748 30 TYPE 110, 

T2796 110 FORMAT (/) 

T2818 37 GO TO 2 

T2826 END 
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PROG SW 1 ONFOR SYMBOL TABLE, PUSH START 

T9999 SIN 
f9989 S I NF 
f9979 COS 
T9969 COSF 
T9959 ATAN 
T99^9 ATANF 
T9939 EXP 
T9929 EXPF 
T9919 LOG 
T9909 LOGF 
T9899 SQRT 
T9889 SQRTF 
T9879 ABS 
T9869 ABSF 
T9859 HEAD 
T9849 HEADF 
T9839 X T9349 
T9339 P T8849 
T8839 OP T8749 
T8739 OY T8649 
T8639 EN T8149 
T8139 EPSLN 
T8129 T0000000DI 
T8119 5100 
T8109 0100 
T8099 A1 
T8089 A2 
T8079 A3 
T8069 A4 
T8059 A5 
T8049 0101 
T8039 0101 
T8029 M 
T8019 0001 



T8009 I 
1799? 



/999 0102 
T7969 0102 
T7979 0002 
T7969 0103 
T7959 0103 
T7949 N 
17939 0003 
T7929 0104 
1791 9 0104 
17909 D 
T7899 R 
T7889 000 
17879 5342944800 
T7869 A 

T7859 5000000001 
17849 B 

T7839 0ooooooo55 

T7829 0004 
17819 J 
T7809 SNW 
T7799 SNWX 
T7789 SNWY 
T7779 SNWXY 
T7769 SNWX 2 
T7759 SNWY2 
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T7?i+9 SCO 5 
T7739 7 
T7729 TE 

T7719 Ooi 
T7709 SO 3 

T7699 5ooooooo5o 
T7689 T25OO00052 

T7679 G"006 
T7669 2 

T7659 198942300*0 
T7649 u 

T7639 ToooooocSi 
T7629 131641890*0 
T7619 5023 

T7609 PW 

T7599 5p4 
T7589 T41 4214051 
T7579 5060 
T7569 5009 
T7559 5008 
T7549 Y 
T7539 W 
T7529 XBAR 
T7519 YBAR 
T7509 SXY 
T7^95 SXX 
T7489 BNXT 
T7479 5007 
17469 8bOZ 
T7^9 CHI 
T7449 »70 
17*39 5071 
T?k29 5105 

®% 

T7399 5? 06 

m89 3107 
|737|5I?7 

T7369 am 



sx 

7339 T96OO00O0t 
129 OU 

.7319 Dt- 
T7309 DE . 

T7299 1302585151 
17289 DUE 
T7279 OLE 
T7269 §080 
T7259 5081 
T7249 5109 

B233 |103 
229 0D13 
T7219 3035 

f7209 5037 
jm m% 

T7189 5ot5 

T7179 5030 
T??65 51 10 
T7159 3110 

SW 1 OFF TO IGNORE SUBROUTINES, PUSH START 
PROCESSING COMPLETE 2 7 



ABSOLUTE VALUE SUBROUTINE - ASS 





DORG 


500G 


05000 










TD 


*fl9, MEM, 


05000 


25 


05019 


00400 




TF 


•M-23, 19989*20 *N,0 


05012 


26 


05035 


1986$ 




TF 


FAC, *-* 


05024 


26 


00060 


00000 




TF 


*+35, 01 


05036 


26* 


05071 


05035 




SM 


*+23, 2, 010 


05048 


12 


05071 


00002 




TF 


FAC-2, *-* 


05060 


26 


00058 


00000 




CF 


FAC -2 


05072 


33 


00058 


00000 




EB 




05084 


42 


00000 


00000 


FAC 


08 


,60 


00060 


00000 






MEM 


08 


,400 


00400 


00000 






NK 


OS 


,7 


00007 


00000 






N 


ss 


,HN-1 


00006 


00000 







DE3B 

If ABS is not the seventh subroutine in the FORTRAN 
•Bbreetine deck, the symbol NN met be redefined appropriately. 
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ALPHABETIC HEADING SUSROUTIKF - iE^D 





DORG 


5000 




05000 










BC2 






05000 


56 . 


05044 


00200 




RCTY 






05012 


34 


00000 


00102 




WATY 


HDNG , 




05024 


39 


05071 


00100 




B 






05036 


49 


05068 


00000 




DORG 


*-3 




05044 








PUNCH 


WACD 


HDNG.,0 




05044 


39 


05071 


00400 


STFLG 


SF 


51 




05056 


32 


00051 


00000 




BB 






05068 


42 


00000 


00000 




DORG 


*-9 




05070 








HONG 


DAC 


40, P 


LOG ED (P) 


LOWER. 


UPPER 










05071 


40 x 


2 






DAC 


41, 


ED (P) 


LOWER 


UPPER 


<? 










05151 


41 x 


2 





DEND 



If the program is recompiled for a paper tape machine, 

the instruction at 05044 must be changed to 39 05071 00200. 



■ O 

LISTINGS 07 EXTRA SUBROUTINES AS INCLUDED IX TORTRAS SUBROUTINE 

DECK 



ABS - 7th Subroutine 



10 j ! 

2505019004001 760503519869 260006000000 
2600058000001 330005800000 42* 
00086 I 



| 10000 

260507105035 T20507 100002 #0010500005060010001 

0010506005086010002 

; 10003 



HEAD - 8th Subroutine 



11 



11000 



460504400200 1 340000000102 390507100100' 4905068000001 f 05000050480 11001 

390507100400 | 320005100000 420000000000 + 0504405080011002 

i ; ' 

000057000000 1005356470045 442457040000 000000005356; 664559000000 !fl 0507005130011003 

000000000064 j 575745590000 000000000045 442457040000 j 000000000000 ; *1 0513005190011004 



535666455900 
00230 



000000000000 645757455900 OOOOf 



| 1 0519005230011005 
I 11006 



20 



